Parity nonconservation in ytterbium ion 
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We consider parity nonconservation (PNC) in singly ionized ytterbium (Yb""") arising from the neu- 
tral current weak interaction. We calculate the PNC electric dipole transition amplitude (EIpnc) 
and the properties associated with it using the relativistic coupled-cluster theory. Elp^c for the 
[4/^*] ^6s — >■ [4/^^] "^5(^3/2 transition in Yb"'" has been evaluated to within an accuracy of 5%. The 
improvement of this result is possible. It therefore appears that this ion is a promising candidate 
for testing the standard model of particle physics. 
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Parity nonconservation (PNC) has been observed in 
a number of neutral atoms; the latest being Yb [l|, [2]- 
The combined experimental and theoretical results for 
Cs PNC is the most accurate to date Q. It is in agree- 
ment with the standard model (SM) of particle physics. 
The experiment on Cs PNC also reports the observa- 
tion of the nuclear anapole moment (NAM) [3, |j] . For 
more stringent tests of the SM, either the accuracy of 
the PNC data from Cs must be improved or very high 
precision measurements should be carried out on other 
candidates. Exploiting the remarkable advances in laser 
cooling and single ion trapping techniques, experiments 
have been proposed to measure PNC in these systems 
[i,.6|. Both Ba+ [i,|i| and Ra+ 0-il are already under 
consideration for such experiments. Moreover, the PNC 
nucleon-nucleon coupling constant related to the NAM 
obtained from Cs PNC and nuclear data do not agree 
[j| , and this clearly calls for further investigation of these 
results. In these circumstances, it would certainly be de- 
sirable to explore PNC in heavy atomic ions Q. Apart 
from Ba+ and Ra+, we had also suggested Pb"*" as a pos- 
sible candidate for PNC measurement [l3] • The study of 
PNC in the [4/"] ^Qs -^ [4/"] 25^3^^ transition in Yb+ 
was suggested by one of us over a decade ago [ll| and it 
is currently being experimentally investigated at the Los 
Alamos National Laboratory [l2{ . 

Assuming that the PNC light-shift in the above tran- 
sition in Yb+ can be measured to sub one per cent ac- 
curacy, it would certainly be desirable to determine the 
accuracy to which ElpNc can be calculated for this tran- 
sition. This knowledge is necessary for choosing a system 
for investigating atomic PNC. In alkali atoms and alkali 



earth-metal ions, the low-lying energy levels which are 
the singly excited states from the ground state usually 
make the dominant contributions to the Elp^c ampli- 
tude. Therefore, a typical sum-over-state approach con- 
sidering only few important singly excited intermediate 
states can give a reasonably accurate result. However, 
Yb"*" has a number of doubly excited low-lying interme- 
diate states, and therefore it would not be prudent to use 
a sum-over-states approach for this ion. 

In this Rapid Communication, we employ the relativis- 
tic coupled-cluster (RCC) method to calculate ElpMc 
using an approach that circumvents summing over inter- 
mediate states [ij|. We also present the results of our 
calculations of the excitation energies and the lifetimes 
for the [4/^^] ^6^1/2,3/2 and the electric quadrupole mo- 
ments of the [4/^''] ^6s — >■ [4/^"*] ^5^3/2 transition, which 
is necessary for evaluating the PNC light-shift associated 
with this transition in Yb+. We analyze the role of the 
electron correlation effects in Elpjyc and compare them 
with those in Ba+ ^ and Ra+ 0. 

In Fig. m we present a diagram of the low-lying en- 
ergy levels of Yb"^ mentioning lifetimes of the important 
states and transitions that can be induced by lasers for 
the measurement of PNC. As shown in the fi gure, t he 
lifetimes of the 5^3/2,5/2 states are fairly large |14l - [l6| . 

The largest contribution to Elp^c in the [4/^"*] ^6s — >■ 
[4/^ ] ^5(i3/2 transition comes from the nuclear spin inde- 
pendent (NSI) neutral weak current interaction between 
the electrons and the nucleus due to the exchange of Zo 
boson. The Hamiltonian for this interaction is given by 
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where Gf{— 2.219 x 10^ in au) is the Fermi constant 
that quantifies the strength of the weak interaction, pnuc 
is the nuclear density which is evaluated using the Fermi 
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FIG. 1: (color online) A schematic diagram of energy levels in 
Yb"*" with transitions shown in red lines that can be induced 
by different lasers for the measurement of parity nonconser- 
vation effect. Lifetimes of different states are denoted by r. 



distribution, 75 is the Dirac matrix, N is the neutron 
number of the atomic system and Qw is known as the 
weak nuclear charge. As the strength of the weak in- 
teraction is very small compared to the electromagnetic 
interactions in an atomic system, it can be treated as a 
first order perturbation. The wave function of an atomic 
state with a single valence electron v can be expressed as 



|*„) = |vI/lo))+Gf !*(!)), (0.2) 

with Gp representing the perturbed parameter, \'i'v ) 
and 1^1) ) are the wave functions corresponding to the 
the Dirac-Coulomb Hamiltonian {Hdc) and its first or- 
der correction due to the PNC weak interaction, respec- 
tively. In our RCC approach, we write 



\^„) - 



e^\l 



S,Mr. 



(0.3) 



where |$^) is the Dirac- Fock (DF) wave function ob- 
tained by appending the valence electron v to the closed- 
shell configuration reference state |$o) that is [4/-'^"'] in 
the present system and T and S^ are the core-virtual and 
valence-virtual excitation operators, respectively. In this 
work, we have considered all possible single and double 
excitations (CCSD method) by taking the operators 



T = Ti +T2 and Sy = 5i„ + S; 
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where subscripts 1 and 2 stand for the single and double 
excitations, respectively. We also consider contributions 
from triple excitations involving the valence electron by 
defining the following operator 
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where e represents the orbital energy and p, q, r and a, b 
represent the virtual and occupied orbitals, respectively. 
When 5*31, is used only to improve the CCSD wave func- 
tions, it is referred to as the CCSD(T) method and when 



its contributions are also estimated then it is known as 
lo-CCSDvT method ^. Further, we separate the T and 
Sv operators to represent the unperturbed and perturbed 
wave functions as 

T = t'^o)^q^j.{i) ^^j Sy^Si"^ +GFSi^\ (0.6) 

where the superscripts and 1 denote unperturbed and 
perturbed operators, respectively. In an explicit form, 
we can write 
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and 1*1,1)) == 



5f)]|<i>.) (0.7) 

^""[T(i)(H-S(°))-K5W]|<i>„). (0.8) 



The RCC equations for the correlation energy and the 
cluster amplitudes equations for both the unperturbed 
and perturbed wave functions have been given in some 
of our earlier papers (for example see [Ij, [l7|). 

An atomic property can be determined in the RCC 
framework by evaluating the matrix element of an oper- 
ator O corresponding to that property 
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Expectation values are 



where Oq 

S'^°^]|$„) with A/'o 
determined by taking f = i. 

The expression for ElpNc between |*/) and |*i) after 
keeping terms up to linear in Gp is given by 
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for the electric dipole (El) operator D, which in the RCC 
approach leads to 
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In TableU we present results for different properties us- 
ing our CCSD(T) method. The results for the appropri- 
ate properties can be used for estimating the uncertainty 
of our Elpj^c result. We also compare our results with 
previous calculations and experimental results wherever 
possible. We analyze the results briefly below. 

Basis functions: We have considered up to 17 s, p, d 
and / and 14 g orbitals using Gaussian type orbital basis 
functions. A Fermi-nuclear charge distribution is used to 
determine the nuclear potential and density. 

Excitation energies: We have given excitation energies 
(EEs) in TableU corresponding to the transitions that are 
important in obtaining an accurate value for Elpj^c re- 
sult and the lifetimes of excited states. Among them EE 



TABLE I: Comparison of different properties in ^^^Yb^ from 
the present and other works (given up to second decimal 
place). Uncertainties are given inside the parentheses. 



TABLE IL ElpNC results in Yb+ in 10""ieao(-Qw'/iV) us- 
ing different approximations. 



Properties 


This work 


Others 


Experiment 


Transitions 


Excitation energies (in cm~^) 


6Si/2 -> 6pi/2 


28109(1000) 


28048" 


27061.82" 


6Sl/2 -^ 6p3/2 


31604(800) 


31411" 


30392.23'' 


6Si/2 -^ 7pi/2 


63518(200) 


63227" 


63706.25' 


6Pl/2 ^ 5^3/2 


4879(800) 


6810" 


4101.02** 


6P3/2 — >■ 5^3/2 


8375(950) 


10173" 


7431.43' 


6P3/2 — >■ 5^5/2 


7015(1000) 


8962" 


6059.54' 


7p3/2 -*■ 5^3/2 


41502(1500) 


43175" 


42633.30' 




El matrix elements (in 


an) 


6Si/2 -^ 6pi/2 


2.72(1) 


2.68/2.73", 2.76^= 




6Si/2 -> 6p3/2 


3.83(1) 


3.77/3.84", 3.87^= 




6Sl/2 -J- 7p3/2 


0.18(2) 


0.15/0.04" 




6Pl/2 -^ 5^3/2 


3.06(2) 


2.97/3.78", 3.19" 




6P3/2 ^ 5^3/2 


1.35(2) 


1.31/1.55", 1.40" 




6p3/2 -> 5^5/2 


4.23(3) 


4.12/4.77", 4.39" 




7Pl/2 -i> 5^3/2 


0.17(2) 


0.08/0.12" 






E2 matrix elements (in 


au) 


6Si/2 — > 5^3/2 


10.2(5) 






6Sl/2 -> 5d5/2 


12.9(7) 




12(1) 


States 


Ahyf constants (in MHz) 


6Sl/2 


13332(1000) 


13172", 12730" 


12645(2)" 


6pi/2 


2516(400) 


2350", 2317^* 


2104.9(1.3)^* 


6p3/2 


322(20) 


311.5", 391^* 


877(20)"= 


7Pl/2 


861(50) 


807.4" 




7p3/2 


123(15) 


110.8" 




5^3/2 


447(20) 


400.48-'^ 


430(43)5 


54/2 


-48(15) 


-12.58^ 


-63.6(7)'' 
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of the 6si/2 —T' Qpi/2 transition contribute significantly 
to the ElpMc result which is obtained to an accuracy 
of around 4% in this work. All our results are compared 
with MBPT(3) results of Safronova and Safronova [1^ 
and also with experimental values [19| . 

El matrix elements: It is important to find the ac- 
curacy of the El matrix elements that contribute to 
ElpNC- We have presented the important El matrix 
elements in the above table. Our results compare rea- 
sonably well with previous calculations fw, i20] . In 18] , 
two different approaches were considered to obtain these 
results which differ from each other substantially. We 
have compared our results with Maui and Angom (their 
published results |20i] have been corrected ^21,]). We com- 
bine our El matrix elements with the experimental wave- 
lengths to determine the lifetimes of the important inter- 
mediate states; 6p ^Pi/2 and 6p ^Pzi2 states and compare 



"^-^^^^^^^^ Method 
Result ^~~~~-~-~~.,_^ 


DF CCSD CCSD(T) lo-CCSDvT 


ElpNC 


6.5311 8.5126 8.4702 8.4704 



20.70 au and 1.15 au, respectively for the 6^3/2 -> ^Pi/2 
transition, we find lifetimes of the 6p ^Pi/2 and 6p ^P^/2 
states as 6.70 ns and 4.75 ns, which agree reasonably 
with the experimental results 7.1(4) ns and 5.5(3) ns of 
Blagoev et al Q, respectively. Also, the branching ratio 
of the decay rate of the 6p ^Pi/2 state to the 5d ^I?3/2 
state has been reported as 0.0483 with 9% uncertainty 
[15| and it agrees well with our result 0.0439. 

E2 matrix elem,ents: We have also determined the E2 
matrix elements of the 6si/2 -^ 5^3/2 and 6si /2 — >■ 5d 



^5/2 



them with the measured values in order to test their accu- 
racies. Using our calculated E2 and Ml amplitudes; i.e. 



transitions which are needed in evaluating the PNC light 
shift. No other results are available for comparison. Tay- 
lor et al have measured the lifetimes of the 5d ^-05/2 state 
as 7.2(3) ms with a branching ratio to the ^-^7/2 state 
of 0.83(3) [16] . From their measurements, we estimate 
the E2 matrix element of the 6si /2 ~> 5^5/2 transition as 
12(1) au and our result is in agreement with that. There- 
fore, we expect that our E2 transition amplitude for the 
6si/2 — > 5(^3/2 transition will also be of similar accuracy. 

Hyperfine structure constants: We have presented the 
magnetic dipole hyperfine structure constants (Ahyf) of 
the low- lying singly excited states of ^^^Yb+. They are 
also compared with other calculations and experimen- 
tal results. Our results are comparatively larger than 
the experimental values, but all the measurements are 
relatively old. A large disagreement between experimen- 
tal and all theoretical results of the 6p ^i-3/2 state calls 
for verification of these results using new experimental 
techniques. However, the accuracies of our results for 
obtaining a reasonable estimate of Elpiqc are currently 
acceptable. For high precision results, it would be neces- 
sary to improve the accuracies of the calculations. 

ElpNC result: We turn here to our Elp^c calculation 
of the 6si/2 — > 5(^3/2 transition where the correlation 
effects seem to be significant. It can be seen from Table 
|TT]that the Elp^c result increases by 30% and 29% at 
the CCSD and CCSD(T) levels, respectively. So there are 
some cancellations between the correlation contributions 
coming via the CCSD and CCSD(T) methods. 

We now proceed to discuss the trends of the correla- 
tion effects by considering the contributions from differ- 
ent RCC terms. In Table Hill we give contributions from 
the DF and individual CCSD(T) terms. By relating the 
important CCSD(T) terms to their corresponding lower 
order MBPT counterparts, we can identify the correla- 
tions effects that are significant. As seen in Table lllli 
three different classes of RCC terms are of crucial impor- 
tance in arriving at the final result. The 'DqT^ term 



TABLE III: Contributions (Contr.) to EIpnc result (in 
xlO-^^ieao{-Qw/N)) in Yb+ at the DF and CCSD(T) 
methods as initial (i — 6s '^Si/2) and final (/ — 5d ^1)3/2) 
perturbed states. Results given as Norm and Others refer to 
contributions from normalization of the wave functions and 
RCC terms those are not mentioned explicitly, respectively. 
Subscripts c and vir represent contributions from core and 
virtual orbitals, respectively. 



Method Initial state 


Contr. 


Final state 


Contr. 


DF 


{HpNcD)c 


0.6441 


{DHpNc)c 


-3.3 X 10"^ 




{DHpNc)vir 


5.8870 


(HpNcD)vir 


-4.3 X 10~^ 


RCC 


rf^'+A, 


0.6745 


OoTi*" 


0.0070 




^o5lP 


7.1973 


S[';'Vo 


1.3424 




^o4^' 


-0.2929 


sifVo 


-0.2022 




slf©o5lP 


0.5115 


s[)^'VoS^^ 


-0.0788 




S^fOo^lP 


-0.3682 


si^'^DoS'^' 


-0.0617 




SiT^oSi]^ 


0.0004 


sifVoS^^' 


-0.0165 




s'fDo^i^) 


0.0271 


si'^'VoS^^' 


0.0096 




Others 


-0.6261 




0.7078 




Norm 


-0.2911 




-0.0670 



represents a sub class of core- valence correlations, whose 
contribution in the present case is slightly larger than DF 
result that is denoted as {DHpnc)c- As we have previ- 
ously discussed [ij,[l3], I^O'S'jJ effectively considers con- 
tributions from singly excited states with {DHpis[c)vir 
as its leading term and involves contributions from im- 
portant pair-correlation and a class of core-polarization 
effects. Clearly this contribution through the perturbed 
5d ^D^/2 state at the DF level is negligible, but after 
including electron correlation effects, its contribution be- 
comes large. The other important contributions com- 
ing via I?0'5'2t, are from the core polarization effects in- 
volving doubly excited states such as 4/^'^ {^F^/2) 6s^ 
and 4/^'^ (^-^7/2) 5(i6s. A sizable amount of contribution 
comes collectively from certain higher order RCC terms 
given as Others. Contributions from the normalization 
of the wave functions (given as Norm) are also not small. 
A detailed comparative study of the trends of the corre- 



lation effects in EIpnc amplitudes of the Si/2 — >• 5^3/2 
transitions in Yb"*" and similar transitions in Ba"*" and 
Ra+ will be reported elsewhere. However, there is an im- 
portant point that we would like to highlight here. There 
is a cancellation between DqSiJ and its complex conju- 
gate (cc) term in Ba"*" and Ra"*", but in contrast in the 
present case it adds up. Therefore, the net correlation 
contribution is larger in Yb"*" than other two ions. 

Analysing results from Table Hill we find that the lead- 
ing contributions from DqT^^' + cc through the core elec- 
trons, [X>o5'jP + sfJ^VaS'^^^] + cc as in the form of the 
single excitations of the valence electrons and the double 
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excitation terms [PoS'^i + S'^fVoS^^' + S)^fVoS^'] + cc 
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are around 8%, 101% and 9% (with opposite sign), re- 
spectively. The uncertainties in the contributions from 
the core and doubly excited states are taken as the dif- 
ferences of the CCSD and lo-CCSD(T) results. Uncer- 
tainties from the singly excited states are estimated from 
the errors of various physical quantities which are given in 
Table U] and taking into account their contributions with 
appropriate weights. With those uncertainties, we esti- 
mate the ElpNc amplitude for the 6s ^S'1/2 -^ 5d ^D^i2 
transition in Yb+ to be 8.5(5) x 10"" iGFeafi{~Qw /N). 
The results can be improved by including the triple and 
quadruple hole-particle excitations in the framework of 
general-order RCC theory [28j . 

In conclusion, Yb"*" appears to be a promising candi- 
date for the study of atomic PNC. By employing the 
RCC method in the singles, doubles and partial triples 
approximation, we have estimated various atomic prop- 
erties including the Elp]\jc amplitude. Our preliminary 
investigation suggests that it is feasible to obtain the- 
oretical results in this system with high accuracy after 
suitable modifications to our present method which has 
yielded a result that is about 5% accurate. A different 
trend for correlation effects is found in the present case 
than Ba+ and Ra"*"; the two ions that were previously 
studied using the same method. 

We thank J. R. Torgerson for communicating with us. 
The computation of the results presented in this work 
were obtained using PRL HPC 3TFL0P cluster. 
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